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Abstract: Grazing is a main human activity in the grasslands of Xinjiang. It is vital to identify the effects of grazing 
on the sustainable utilization of local grasslands. However, the effects of grazing on net primary productivity (NPP), 
evapottanspiration (ET) and water use efficiency (WUE) in this region remain unclear. Using the spatial Biome- 
BGC grazing model, we explored the effects of grazing on NPP, ET and WUE across the different regions and 
grassland types in Xinjiang during 1979-2012. NPP, ET and WUE under the grazed scenario were generally lower 
than those under the ungtazed scenario, and the differences showed increasing trends over time. The decreases in 
NPP, ET and WUE varied significantly among the regions and grassland types. NPP decreased as follows: among 
the regions, Northern Xinjiang (16.60 g C/(m*a)), Tianshan Mountains (15.94 g C/(ma)) and Southern Xinjiang 
(-3.54 g C/(m2a)); and among the grassland types, typical grasslands (25.70 g C/(m?-a)), swamp meadows (25.26 g 
C/(m2a)), mid-mountain meadows (23.39 g C/(m2a)), alpine meadows (6.33 g C/(m2a)), desert grasslands (5.82 g 
C/(m?-a)) and saline meadows (2.90 g C/(m2-a)). ET decreased as follows: among the regions, Tianshan Mountains 
(28.95 mm/a), Northern Xinjiang (8.11 mm/a) and Southern Xinjiang (7.57 mm/a); and among the grassland types, 
mid-mountain meadows (29.30 mm/a), swamp meadows (25.07 mm/a), typical grasslands (24.56 mm/a), alpine 
meadows (20.69 mm/a), desert grasslands (11.06 mm/a) and saline meadows (3.44 mm/a). WUE decreased as 
follows: among the regions, Northern Xinjiang (0.053 g C/kg H20), Tianshan Mountains (0.034 g C/kg H20) and 
Southern Xinjiang (0.012 g C/kg H20); and among the grassland types, typical grasslands (0.0609 g C/kg H20), 
swamp meadows (0.0548 g C/kg H20), mid-mountain meadows (0.0501 g C/kg H20), desert grasslands (0.0172 g 
C/kg H20), alpine meadows (0.0121 g C/kg H20) and saline meadows (0.0067 g C/kg H20). In general, the 
decreases in NPP and WUE were more significant in the regions with relatively high levels of vegetation growth 
because of the high grazing intensity in these regions. The decreases in ET were significant in mountainous areas 
due to the terrain and high grazing intensity. 


Keywords: grazing effect; grassland type; net primary productivity; evapotranspiration; water use efficiency; Biome- 
BGC grazing model 


1 Introduction 


Grazing is common in grasslands (Zuo et al., 2012; Eldridge et al., 2016) and has significant effects 
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on the ecological functions of grassland ecosystems (Dong et al., 2004a; Teague et al., 2011; Zhang 
et al., 2015; Wang et al., 2016a; Xie et al., 2016). The interactions between grassland ecosystems 
and herbivores are highly complex, and many factors are involved. Therefore, grazing has been a 
popular but challenging topic in the field of grassland ecology. However, researches regarding the 
effects of grazing on net primary productivity (NPP), evapotranspiration (ET) and water use 
efficiency (WUE) in different types of grasslands are limit at the regional scale, particularly in arid 
areas. 

NPP, ET and WUE are the three key factors that reflect the functions of grassland ecosystems. 
NPP reflects the production capacity of the grassland ecosystems and is highly sensitive to climate 
change and human activities (Han et al., 2014). ET is the main form of water loss in arid grasslands, 
and the balance between precipitation and ET is crucial for determining ecosystem structure, 
function, and productivity (Eldridge et al., 2016; Liu et al., 2016). WUE reflects how efficiently 
plants use water to produce biomass and is used as an indicator of the suitability of an environment 
for plant growth (Dong et al., 2004b; Niu et al., 2011; Zhang et al., 2016). Furthermore, compared 
with plants in humid grasslands, plants in arid grasslands exhibit unique ecological processes in 
response to water stress (Luo et al., 2012; Liu et al., 2016). NPP, ET and WUE are highly correlated 
with water availability. Therefore, studies on NPP, ET and WUE are particularly important in dry 
areas. 

Grazing can affect plant growth through the removal of a certain amount of vegetation and 
through the deposition of feces on the soil surface (Olejniczak, 2011; Irisarri et al., 2016). Grazing 
can also affect ET through the removal of vegetation (Wang et al., 2016b), and WUE may be 
changed by grazing. Recent studies show that the responses of NPP and ET to grazing differ under 
different environmental conditions (Luo et al., 2012; Zhang et al., 2015; Eldridge et al., 2016; Wang 
et al., 2016b). Therefore, NPP, ET and WUE might show diverse effects across different regions 
and grassland types. Understanding the effects of grazing on NPP, ET and WUE over large areas is 
helpful for making appropriate grazing management decisions in accordance with different 
grassland ecosystems. 

Xinjiang is located in an arid region of the Eurasian hinterland that encompasses widely 
distributed grassland ecosystems, and grazing is the main anthropogenic disturbance in local 
grasslands. With increasing grazing intensity in recent decades, grassland degradation has become 
widespread (Han et al., 2014, 2016; Liu et al., 2016). Therefore, the research on the effects of 
grazing on NPP, ET and WUE is vital for the sustainable utilization of grassland resources in 
Xinjiang, where distinct differences in elevation, climate, grassland type and grazing management 
lead to spatial heterogeneity in NPP, ET and WUE. Therefore, this area is conducive to conducting 
research on the effects of grazing under different environmental and grazing conditions in arid 
areas. 

At present, model simulation is considered the most effective method for studying the effects of 
grazing over large areas (Laniak et al., 2013; Dugdale et al., 2016; Guo et al., 2016) as different 
scenarios can be evaluated using model estimations. Therefore, we can effectively evaluate the 
effects of grazing using model simulations to identify the differences between ‘grazed’ and 
‘ungrazed' scenarios and to compare NPP, ET and WUE among studies conducted in the same area. 
The process-based Biome-BGC grazing model can effectively estimate the effects of grazing on 
the fluxes and storage of carbon and water in grassland ecosystems (Wang et al., 2005; Luo et al., 
2012; Han et al., 2016). 

The objective of this study was to determine the effects of grazing on NPP, ET and WUE in the 
grasslands of Xinjiang during 1979-2012 across different regions and grassland types using the 
spatial Biome-BGC grazing model. 


2 Materials and methods 


2.1 Study area 


Xinjiang, covering an area of 1.66x10° km?, is an important region for the production of livestock. 
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The climate is typically temperate with arid characteristics. The annual mean temperature is 10.4°C 
and the annual mean precipitation is 188 mm. The Tianshan Mountains divide Xinjiang into two 
sub-regions: Northern Xinjiang and Sourthern Xinjiang (Fig. 1a). The mean annual precipitation in 
Southern Xinjiang ranges from 20 to 100 mm, and the annual mean temperature ranges from 10.0°C 
to 13.8°C. In contrast, the mean annual precipitation in Northern Xinjiang ranges from 100 to 500 
mm, and the annual mean temperature ranges from 4.0°C to 8.8°C (Gu et al., 2010; Du et al., 2015; 
Jiapaer et al., 2015). Grasslands account for 34.4% of the territory in Xinjiang and include alpine 
meadows, mid-mountain meadows, typical grasslands, desert grasslands, saline meadows, and 
swamp meadows (Huang et al., 2017), with surface areas of 1.09x10°, 3.14x10°, 3.35x10°, 
3.63x10°, 3.01x10°, and 5.4x104 km?, respectively (Fig. 1b). Most grasslands in Xinjiang are used 
for grazing, but grazing intensity differs significantly among the different types of grasslands (Figs. 
lc and d). With the increase in grazing intensity in recent decades, large areas of grasslands have 
suffered severe deterioration (Xu et al., 2012; Rong et al., 2014). 
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Fig. 1 Distributions of elevation (a), grassland type (b), pasture type (c) and average grazing intensity (d) in 
Xinjiang during 1979-2012. NPP, net primary productivity; ET, evapotranspiration. 


2.2 Data 


The data in the study included model validation data, grazing data, meteorological data, and other 
auxiliary data. Grazing data, meteorological data, and other auxiliary data served as the model input 
data. 


2.2.1 Model validation data 


Model validation data, including NPP and ET, were collected from both field observations and 
previous publications by Li et al. (2013) and Han et al. (2014). There was a total of 48 plots with 
annual NPP data, among which 6 plots were located inside enclosures (i.e., not grazed) and 42 were 
located outside enclosures and were grazed. Because it was difficult to measure NPP intake by 
livestock, all the NPP data collected from the grazed plots represented the NPP remaining in the 
ecosystems after grazing. Daily ET data during 2013-2015 were obtained from field observations 
in a typical grassland in Urumqi County (43°33’N, 87°12’E; 1648 m a.s.l.) using an HL20 Bowen 
ratio system at the center of a 100 mx100 m enclosure without grazing. Daily ET data from October 
2012 to September 2013 were obtained from field observations in an alpine meadow in Aksu 
(41°42’N, 80°10'E; 3550 m a.s.l.), using a microlysimeter (Guo et al., 2015) inside an enclosure 
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without grazing. Daily ET values from 22 July to 2 September, 2010, were collected using a Bowen 
ratio system outside an enclosure in a desert grassland in Qitai County (44°11'N, 89°26’E; 617 m 
a.s.l.) and were also obtained from a previous publication by Yan et al. (2015). 


2.2.2 Grazing data 


Elevation, precipitation and the distance from residential areas were used to distinguish different 
pasture types. Summer pastures were located in mountainous areas above 2500 m and in plain areas 
with a mean annual precipitation of 630 mm. Spring-autumn pastures were located in mountainous 
areas below 2500 m and in plain areas with a mean annual precipitation of 126 mm. Winter pastures 
were typically close to mountainous permanent settlements and in plain areas with light snowfall 
and a mean annual precipitation of 210 mm (Banks and Doman, 2001; Kerven et al., 2011; Han et 
al., 2014). Grazing intensity data at the regional scale, which were generated using interpolation in 
accordance with sampling points of global livestock numbers and environmental data, were 
provided by the Food and Agriculture Organization of the United Nations (FAO) 
(http://www.fao.org/docrep/010/a1259e/a1259e00.htm). We corrected the grazing intensity data in 
accordance with statistics from the Xinjiang yearbook, ensuring a high degree of data precision 
(Jin, 2012). Grazing data at model validation points were collected through field investigations. 
The livestock included cattle, sheep, goats, horses, camels, and yaks. All livestock were converted 
into sheep units in accordance with the standards provided by the Ministry of Agriculture of the 
People's Republic of China (http://www.chinaforage.com/standard/zaixuliang.htm). According to 
these standards, one goat equals 0.9 sheep, one cattle equals 6.0 sheep, one horse equals 6.0 sheep, 
one yak equals 4.5 sheep, and one camel equals 8.0 sheep. Figure 2 shows the inter-annual variation 
of grazing intensity in the Xinjiang grasslands during 1979-2012. 
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Fig.2 Interannual variation in grazing intensity in the grasslands of Xinjiang during 1979-2012 


2.2.3 Meteorological data 


Different meteorological data were required for this study. These include: the daily maximum and 
minimum air temperatures, daily average air temperature, humidity, incident solar radiation, 
precipitation and day length. All these regional data were extracted from the China Meteorological 
Forcing Dataset using C++. The precision of these data was better than that of the existing grid- 
based meteorological data for this region. With a spatial resolution of 10 km, they can be used for 
different terrestrial modeling purposes (http://www.tpedatabase.cn/portal/MetaDataInfo.jsp? 
MetaDatald=249369). The meteorological data at the model validation points were collected 
through field observations. 

2.2.4 Other auxiliary data 


Elevation data were derived from the WorldClim database (Hijmans et al., 2005) 
(http://www.worldclim.org/download). Soil data, including effective soil depth and sand/silt/clay 
percentages, were derived from the harmonized world soil database (HWSD; 
FAO/IIASA/ISRIC/ISS-CAS/JRC, 2012). A map of grassland types in Xinjiang was generated 
based on a map of vegetation types in the area according to elevation and vegetation characteristics. 
Elevation and soil data at the model validation points were obtained through field investigations. 


HUANG Xiaotao et al.: Effects of grazing on net primary productivity, evapotranspiration and water... 


2.3 The Biome-BGC grazing model 


The Biome-BGC grazing model, a process-based ecosystem model, can effectively simulate 
grazing effects on the fluxes and storage of carbon and water in grassland ecosystems over large 
areas. The grazing process in the model is suitable for arid and semi-arid grasslands (Seligman et 
al., 1992; Luo et al., 2012; Han et al., 2016). 

In the model, NPP is estimated by Equation 1 (Luo et al., 2012). 


NPP=C' eg +Giitter +D, » (1) 


litter 


where C'veg (g C/(hm?-d)) is vegetative carbon; Chiter is litter carbon (g C/(hm?.d)); and D., (g 
C/(hm?.d)) is the carbon consumed by animals and is calculated by Equation 2 (Luo et al., 2012). 


D, =G,S, CF z (Ce Jy ) (0<D;<SDx), (2) 


where Ge is the grazing efficiency of the livestock (hm?/(sheep-d)), S; is the grazing intensity 
(sheep/hm?), Ciear is the C in the leaf biomass (g C/m7), (Ciear)u is the residual aboveground Cieaf 
that is unavailable to livestock (g C/m?), and Dx is the consumption rate of the livestock based on 
satiation (g C/(sheep-d)). 

In the model, ET is calculated using the Penman-Monteith equation, as shown in Equation 3 
(Running and Coughlan, 1988; Zhao and Zhao, 2014). 


A(R, -G)+pC, (e, =e, )/r, 
A+y(1+r)/r, 


ENZ ; (3) 
where ET is the actual ET (mm/d), A is the slope of the curve of vapor pressure versus temperature 
(kPa/°C), Rn is the net radiation (MJ/(m2.d)), G is the soil heat flux (MJ/(m2.d)), p is the air density 
(kg/m), Cp is the specific heat of air (MJ/(kg-°C)), es is the saturation vapor pressure (kPa), ea is 
the actual vapor pressure (kPa), y is the psychrometric constant (kPa/°C), rs is the bulk surface 
resistance (s/m), and ra is the aerodynamic resistance (s/m). In the Penman-Monteith equation, 
energy is divided into two components. One component is used in the transpiration of water, and 
the other component is used in evaporation. 
In this study, grassland WUE was calculated by Equation 4 (Han et al., 2013). 


WUE=NPPEIET, (4) 


where WUE is expressed as g C/kg H20. 
Two scenarios, grazing and ungrazed, were considered to evaluate grazing effects on NPP, ET 
and WUE for comparison with studies conducted in the Xinjiang grasslands from 1979-2012. 


3 Results 


3.1 Model validation 


In this study, Biome-BGC grazing model was validated with respect to NPP and ET estimates using 
field data collected through observations and the data from previous publications in Xinjiang 
grasslands. Figure 3 shows that the model accurately estimated the annual NPP, with the simulated 
data corresponding closely to field data from previous publications and the observations under the 
both grazing (R7=0.88, RMSE=26.92) and ungrazed scenarios (R?=0.83, RMSE=36.07). Due to 
difficulties in measuring NPP intake by livestock, the NPP values associated with grazing, 
including observations and simulations, represent the NPP remaining in the ecosystems after 
grazing. 

We conducted the ET simulation validation through comparison with observed ET data from a 
typical grassland in Urumqi County under the ungrazed scenario (R?=0.76, RMSE=0.67), in an 
alpine meadow in Aksu under the ungrazed scenario (R7=0.66, RMSE=0.76), and in a desert 
grassland in Qitai County under the grazing scenario (R?=0.66, RMSE=0.30). The results showed 
that the Biome-BGC grazing model performed well in estimating ET (Figs. 4a—c). 
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Fig. 3 Relationship between simulated and observed annual NPP in the grasslands of Xinjiang under the grazing 
scenario (a) and ungrazed scenario (b). NPP values associated with grazing represent the NPP remained in the 
ecosystems after grazing. 
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Fig. 4 Comparisons of daily ET based on simulated and observed field data for a typical grassland in Urumqi 
County under the ungrazed scenario (a), an alpine meadow in Aksu under the ungrazed scenario (b), and a desert 
grassland in Qitai County under the grazing scenario (c) 


3.2 Impacts of livestock grazing on NPP 


Figure 5 shows the spatial distribution of the differences in the average annual NPP between the 
grazing and ungrazed scenarios in the grasslands of Xinjiang during 1979-2012 and the variation 
in NPP over time. In contrast to NPP under the ungrazed scenario, NPP under the grazing scenario 
decreased in most regions, indicating widespread grassland degradation during 1979-2012. 
Significant decreases in NPP were mainly distributed in the Altay Mountains and part of the 
Tianshan Mountains. Increases in NPP were mainly distributed at middle to low elevations on the 
northern slope of the Kunlun Mountains and in the western part of Southern Xinjiang (Fig. 5a). In 
general, the NPP in the grasslands of Xinjiang during 1979-2012 was lower under the grazing 
scenario than under the ungrazed scenario. In addition, similar inter-annual fluctuations were 
observed between the grazing and ungrazed scenarios. A fluctuating increase in NPP was observed 
under the ungrazed scenario, while a fluctuating decline was observed under the grazing scenario 
(Fig. 5b). 
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Fig.5 NPP simulation results under the grazing and ungrazed scenarios in Xinjiang grasslands during 1979-2012: 
(a) spatial distribution of the average annual NPP difference and (b) interannual variation in NPP. NPP 
difference=NPP grazea—NPPungrazea 


In general, the decrease in NPP under the grazing scenario was significant in Northern Xinjiang 
and the Tianshan Mountains. In contrast, NPP generally increased under the grazing scenario in the 
Southern Xinjiang (Table 1). 

Table 1 Net primary productivity (NPP) statistics for different regions 


IND nera NPP grazed Difference Percentage 
(g C/(m°.a)) (%) 
Northern Xinjiang 192.51 175.91 —16.60 8.62 
Southern Xinjiang 84.48 88.02 3.54 4.19 
Tianshan Mountains 287.25 271.31 —15.94 5.55 


Note: NPPungrazea, NPP under the ungrazed scenario; NPPgrazea, NPP under the grazing scenario. 


The decrease in NPP under the grazing scenario declined in the following order: typical 
grasslands, swamp meadows, mid-mountain meadows, alpine meadows, desert grasslands, and 
saline meadows. In general, the decreases in NPP were more significant in grasslands with 
relatively good vegetation growth, such as swamp meadows, mid-mountain meadows, and typical 
grasslands (Table 2). 


Table 2 NPP statistics for different grasslands 


NPPungrazea NPP grazed Difference Percentage 

(g C/(m?.a)) (%) 
Alpine meadows 148.32 141.99 —6.33 4.27 
Mid-mountain meadows 366.08 342.69 -23.39 6.39 
Typical grasslands 323.05 29755 -25.70 7.96 
Desert grasslands 159.99 154.17 -5.82 3.64 
Saline meadows 81.93 79.03 -2.90 3.54 
Swamp meadows 326.52 301.26 -25.26 7.74 


3.3 Impacts of livestock grazing on ET 


Figure 6 illustrates the spatial distribution of the differences in average annual ET between the 
grazing and ungrazed scenarios in the Xinjiang grasslands over the period 1979-2012 and the 
variation in ET over time. In contrast to ET under the ungrazed scenario, ET under the grazing 
scenario decreased in most regions of Xinjiang. Significant decreases in ET under the grazing 
scenario were mainly observed in part of the Tianshan Mountains. Increases in ET were mainly 
observed in part of Southern Xinjiang and the north part of Northern Xinjiang (Fig. 6a). In general, 
ET was lower under the grazing scenario than under the ungrazed scenario from 1979-2012. In 
addition, inter-annual fluctuations were similar between the grazing and ungrazed scenarios. An 
obvious fluctuating increase in ET was observed under the ungrazed scenario, whereas no 
significant fluctuating increase was observed under the grazing scenario (Fig. 6b). 
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Fig. 6 ET simulation results under the grazing and ungrazed scenarios in Xinjiang grasslands from 1979-2012: 
(a) spatial distribution of average annual differences in ET and (b) interannual variations in ET. ET 
difference=ET grazed—ET ungrazed. 


The decrease in ET under the grazing scenario was the most significant in the Tianshan 
Mountains. Decreases in ET were slightly larger in Northern Xinjiang grasslands than in Southern 
Xinjiang grasslands. However, due to the lower ET values in the grasslands of Southern Xinjiang, 
the percentages by which ET decreased were slightly higher in Southern Xinjiang than in Northern 
Xinjiang (Table 3). 

Table 3 ET statistics for different regions 


iima B Difference Percentage 
(mm/a) (%) 
Northern Xinjiang 221.31 213.20 —8.11 3.67 
Southern Xinjiang 183.22 175.65 -7.57 4.13 
Tianshan Mountains 357.94 328.99 -28.95 8.08 


Notes: ETungrazea, ET under the ungrazed scenario; ET grazea, ET under the grazing scenario. 


The decreases in ET declined in the following order: mid-mountain meadows, swamp meadows, 
typical grasslands, alpine meadows, desert grasslands and saline meadows. The percentages by 
which ET decreased declined in the following order: alpine meadows, mid-mountain meadows, 
typical grasslands, swamp meadows, desert grasslands and saline meadows. In general, the 
decreases in ET were more significant in mountainous grasslands, such as alpine meadows, mid- 
mountain meadows and typical grasslands (Table 4). 


Table 4 ET statistics for different grasslands 


imeri Bi Difference Percentage 
(mm/a) (%) 
Alpine meadows 229.70 209.04 -20.69 9.01 
Mid-mountain meadows 387.70 358.37 -29.30 7.56 
Typical grasslands 343.90 319.33 -24.56 7.14 
Desert grasslands 236.20 225.06 —11.06 4.72 
Saline meadows 154.70 158.15 3.44 2p 
Swamp meadows 358.30 333.26 -25.07 6.70 


3.4 Impacts of livestock grazing on WUE 


Figure 7 illustrates the spatial distribution of the differences in average annual WUE between the 
grazing and ungrazed scenarios in the grasslands of Xinjiang over the period 1979-2012 and the 
variation in WUE over time. In contrast to the WUE under the ungrazed scenario, the WUE under 
the grazing scenario decreased in most Xinjiang grasslands. And significant decreases in WUE 
were mainly distributed in the northwestern Tianshan Mountains and the eastern Altay Mountains. 
WUE increased mainly in the Altun Mountains, the southwestern Tianshan Mountains, and the 
northern slope of the middle and low elevations of the Kunlun Mountains (Fig. 7a). The inter- 
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annual fluctuations in WUE were lower under the grazing scenario than under the ungrazed 
scenario during 1979-2012. In addition, an obvious fluctuating decline in WUE was observed 
under the grazing scenario in Xinjiang, whereas under the ungrazed scenario no significant 
fluctuating decline was observed (Fig. 7b). 
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Fig. 7 Water use efficiency (WUE) simulation results under the grazing and ungrazed scenarios in Xinjiang 
grasslands from 1979-2012: (a) spatial distribution of average annual WUE difference and (b) interannual variation 
in WUE. WUE difference=WU Egrazea~W U Eungrazed. 


Under the grazing scenario, the decrease in WUE was most significant in the Northern Xinjiang 
grasslands and was least significant in the Southern Xinjiang grasslands. Decreases in WUE in the 
Tianshan Mountain grasslands were also obviously higher than those in the Southern Xinjiang 
grasslands. However, due to the higher WUE in the Tianshan Mountains, the percentages by which 
WUE decreased were slightly higher in the Tianshan Mountains than in the Southern Xinjiang 
grasslands (Table 5). 


Table 5 Water use efficiency (WUE) statistics for different regions 


WUE ngrazea WUE grazea Difference Percentage 
(g C/kg H20) (%) 
Northern Xinjiang 0.730 0.677 —0.053 7.26 
Southern Xinjiang 0.260 0.248 —0.012 4.62 
Tianshan Mountains 0.690 0.656 —0.034 4.88 


Note: WUEungrazea, WUE under the ungrazed scenario; WUEgrazea, WUE under the grazing scenario. 


WUE declined in the following order: typical grasslands, swamp meadows, mid-mountain 
meadows, desert grasslands, alpine meadows and saline meadows. In general, the decrease in WUE 
was more significant in grasslands with relatively high levels of vegetation growth, such as swamp 
meadows, mid-mountain meadows and typical grasslands (Table 6). 


Table 6 WUE statistics for different grasslands 


WUE ngrazea WUE grazea Difference Percentage 
(g C/kg H20) (%) 
Alpine meadows 0.3700 0.3579 0.0121 A27 
Mid-mountain meadows 0.8600 0.8099 —0.0501 5.83 
Typical grasslands 0.8400 0.7791 —0.0609 IA 
Desert grasslands 0.3400 0.3228 —0.0172 5.06 
Saline meadows 0.5200 0.5133 —0.0067 129) 
Swamp meadows 0.7700 0.7152 —0.0548 Til 


4 Discussion 


4.1 Uncertainty analysis 


It was confirmed that the grazing process modeling with Biome-BGC grazing is suitable for 
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grasslands in arid and semi-arid areas (Seligman et al., 1992; Luo et al., 2012). This model can 
effectively estimate the carbon emissions of grasslands in Central Asia (Han et al., 2016). In this 
study, we conducted model validation with respect to NPP and ET simulation results in the 
grasslands of Xinjiang to confirm the reliability of the model. However, the grazing process is 
complex, and some components were not fully considered in the model. As with many other 
models, uncertainty in the results of the model used here is inevitable. For example, the effects of 
trampling remain unclear and are not effectively considered in the Biome-BGC grazing model, 
which inevitably leads to uncertainty in the results. Ecological processes can be influenced in 
grassland ecosystems due to the soil density change by trampling (Neff et al., 2005; Bell et al., 
2011). We hope that the grazing process of the model can be improved in the future to effectively 
consider trampling effects. 

In addition to the uncertainty associated with the model itself, uncertainty in the model results 
can also arise from the input data. Meteorological data are the major drivers that influence the 
model results (Bond-Lamberty et al., 2005; Jung et al., 2007). The meteorological data used in this 
study were extracted from the China Meteorological Forcing Dataset 
(http://www.tpedatabase.cn/portal/MetaDataInfo. jsp? MetaDatald=249369). The precision of these 
data was higher than that of other existing meteorological grid datasets for this region. Grazing data 
were the key inputs used to simulate the effects of grazing (http://www.fao.org/docrep/010/a1259e 
/a1259e00.htm) and were provided by FAO in this study. We corrected the data in accordance with 
Statistics from the statistical yearbook, which ensured high data precision (Jin, 2012). 


4.2 Effects of livestock grazing on NPP, ET and WUE in grasslands 


The results suggest that NPP decreased as a result of livestock grazing in the grasslands of Xinjiang 
during 1979-2012, indicating that unsustainable grazing, as a result of increased grazing intensity, 
and widespread grassland degradation has occurred in Xinjiang in recent decades (Mikola et al., 
2009; Su et al., 2017; Wang et al., 2017). However, the increases in NPP in some regions indicated 
that moderate grazing may promote grass growth (Paige, 1992; Nilsson et al., 1996; Chen et al., 
2007; Rótolo et al., 2007). The effects of grazing on NPP varied among the different regions and 
grassland types due to differences in the climate and grazing intensity. In general, decreases in NPP 
were more significant in the grasslands with high levels of vegetation growth, which can mainly be 
attributed to the high grazing intensity in these grasslands (Fig. 1). 

Our study shows that the difference in ET between the grazing and ungrazed scenarios increased 
due to the increased grazing intensity in recent decades. In plain grasslands, the decrease in ET was 
mainly caused by decreased transpiration. Grazing can remove parts of leaf blades and decrease 
transpiration accordingly (Verón et al., 2011; Eldridge et al., 2016). The decrease in ET in 
mountainous grasslands was more significant than the decrease in ET in plain grasslands. In 
general, the grazing intensity was higher in mountainous grasslands than in plain grasslands, and a 
higher grazing intensity may lead to lower transpiration. Grazing led to the destruction of 
vegetation, which not only decreased transpiration but also generated more runoff (Leitinger et al., 
2010; Orr et al., 2016). Therefore, decreases in ET were more significant under the grazing scenario 
in the mountainous grasslands than in the plain grasslands. 

WUE under the grazing scenario was lower than that under the ungrazed scenario. Furthermore, 
a strong fluctuating decline in WUE was observed under the grazing scenario, whereas no such 
pattern was observed under the ungrazed scenario. Therefore, increased grazing intensity led to a 
remarkable decrease in WUE. Overgrazing can remove excessive amounts of leaves (the main 
photosynthetic organ of plants) and is not conducive to grass growth (Zhao et al., 2007; Eldridge 
et al., 2016). In recent decades, overgrazing has been widespread and has caused decreases in WUE 
in most regions of Xinjiang. In general, decreases in WUE were more significant in the grasslands 
with high levels of vegetation growth because these regions were generally associated with high 
grazing intensities. 


4.3 Comparisons with previous studies 


The effects of grazing on the grasslands of Xinjiang have been studied using a variety of methods. 
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Using observational data, Long et al. (2010) conducted researches on evaporation from grassland 
soil under different rainfall intensities in Qitai County of Xinjiang. Yan et al. (2013) explored the 
relationships of ET with multiple environmental factors in an Achnatherum splendens grassland in 
Xinjiang during summer and autumn. Using a Biome-BGC grazing model, Luo et al. (2012) 
estimated aboveground net primary productivity (ANPP) at four sites in the northern Tianshan 
Mountains under different grazing intensities. This was the first use of this grazing process model 
to investigate the effects of grazing in the Xinjiang grasslands. The authors disclosed that ANPP 
can be promoted under moderate grazing scenario for grasslands under water stress. However, the 
previous studies which based on sites did not clearly provide an overall understanding of these 
phenomena over large scale areas, particularly in areas with complex surfaces, such as the 
grasslands of Xinjiang. Han et al. (2014) conducted a study on carbon emission from grasslands 
using the spatial Biome-BGC grazing model and presented the first carbon budget estimate at the 
regional scale for Xinjiang. However, the grazing intensity in that study was obviously lower than 
the real grazing intensity, leading to bad accuracy in terms of results. In this study, to ensure a high 
precision of the grazing data, we adjusted the grazing intensity data in accordance with information 
from the Xinjiang statistical yearbook (Jin, 2012). In addition, the effects of grazing on NPP, ET 
and WUE in Xinjiang were estimated for different regions and grassland types using the spatial 
Biome-BGC grazing model to investigate the spatial differences in grazing effects in detail. 


5 Conclusions 


NPP, ET and WUE of the Xinjiang grasslands under the grazing scenario were generally lower than 
those under the ungrazed scenario during 1979-2012. In addition, the differences between grazing 
and ungrazed scenarios increased as a result of increased grazing intensity over time. Grazing 
turned the NPP from the increased volatility to the decreased volatility in the grasslands of Xinjiang 
between 1979 and 2012. Compared with the ungrazed scenario, grazing led to a slower increase in 
ET and a more significant decrease in WUE over the study period. 

Strong regional variations occurred in NPP, ET and WUE differences between the grazing and 
ungrazed scenarios among the regions and grassland types due to variations in climate and grazing 
management. Widespread overgrazing has contributed to decreases in NPP in most regions in recent 
decades. Because of the decrease in transpiration caused by grazing, ET decreased in the plains 
regions. In mountainous grasslands, in addition to a decrease in transpiration, an increase in runoff 
caused by grazing was an important contributor to the decrease in ET. WUE decreased in most 
regions in conjunction with the destruction of vegetation caused by overgrazing. This study 
suggests that grazing had a significant impact on carbon and water cycling in Xinjiang grassland 
ecosystems during 1979-2012. 
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